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CL|| Theoretical results for the cross section of heavy quark production near the threshold at NNLO of 
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> 

X 



Talk given at 10th Lomonosov conference on Elementary Particle Physics, 
August 2001, Moscow, Russia 



1 Introduction 



The heavy quark production near thresholds will be thoroughly investigated at future accelerators, 
e.g. H]. One can study the qq systems near threshold in e^e~ annihilation H, |^ and 77 collisions [Q, |5|. 
In e^e~ annihilation the production vertex is local (the electromagnetic current in case of photon 
and/or neutral weak current in case of Z-boson), the basic observable is a production cross section 
which is saturated by S-wave (for the vector current). In 77 collisions the production vertex is nonlocal 
and both S- and P-waves can be studied for different helicity photons, the number of observables is 
larger (cross sections as, crp, S-P interference). In e~^e~ annihilation high precision data are already 
available for bb and expected for tt. The 77 — > qq experiments are planned. 



2 Features of theoretical description 

Heavy quarks near the production threshold moves slowly that justifies the use of the nonrelativistic 
quantum mechanics for their description [6-8] . Being much simpler than the comprehensive relativistic 
treatment of the bound state problem with Bethe-Salpeter amplitude Q, the nonrelativistic approach 
allows one to take into account exactly such an essential feature of the near-threshold dynamics as 
Coulomb interaction For unstable heavy quarks with a large decay width it is possible to compute 
the cross section near threshold point-wise in energy because the large decay width suppresses the long 
distance effects of strong interaction The ti-pair near the production threshold is just a system 
that satisfies the requirement of being nonrelativistic. Therefore the description of tt-system near the 
production threshold ^/s ~ 2mt (y/s is a total energy of the pair, mt is the top quark mass) is quite 
precise within nonrelativistic QCD (NRQCD). Reasons for this accuracy are related to the large mass 
of the top: 



The top quark is very heavy rrit = 175 GeV [12| and there is an energy region of about 8 — 10 GeV 
near the threshold where the nonrelativistic approximation for the kinematics is very precise. 
For ^/s = 2mt + E with < 5 GeV the quark velocity is small 



Relativistic effects are small and can be taken into account perturbatively in v (even in u^). 

• The strong coupling constant at the high energy scale is small as{mt) ~ 0.1 that makes the per- 
turbative mapping of QCD onto the low energy effective theory (NRQCD) numerically precise. 

• The decay width of top quark is large, Tt = 1.43 GeV; infrared (small momenta) region is 
suppressed and PT calculation for the cross section near the threshold is reliable point-wise in 
energy. 

Because ~ v and the ratio as/v is not small, the Coulomb interaction is enhanced. The ordinary 
perturbation theory for the cross section (with free quarks as the lowest order approximation) breaks 
down and all terms of the order (as/v)'^ should be summed. The expansion for a generic observable 
f{E) in this kinematical region has the form 

f{E) = Uas/v) + ash{as/v) + alf2{as/v) + . . . (2) 
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where fi{as/v) are some (not polynomial) functions of the ratio ag/v, fo{ois/v) is a result of the pure 
Coulomb approximation (or a kind of its improvement). The expansion in Og in eq. (|2|) takes into 
account the perturbative QCD corrections to the parameters of NRQCD and relativistic corrections 
(in the regime v ~ a^). 

For the e^e~ tt process mediated by the photon the NNLO analysis is well known. The basic 
quantity is the vacuum polarization function 

n(^) = i I {Tjemix)jemme'i''dx, = {2mt + E f . (3) 

Near the threshold (for small energy E) NRQCD is used. The cross section is saturated with S-wave 
scattering. In this approximation the polarization function near the threshold to the NNLO accuracy 
in NRQCD is given by 

27r 

U{E) = —Ch{ag)Co{E/mt)G{E- 0, 0) . (4) 



The pole mass definition is used for rrit (e.g. ||l^), ag is the strong coupling constant. Ch{o-s) is the 
high energy coefficient. G{E] 0, 0) is the nonrelativistic Green function (GF). The quantity Co{E /mt) 
describes the contributions of higher dimension operators within the effective theory approach. These 
contributions have, in general, a different structure than the leading term. To the NNLO of NRQCD 
the contribution of higher dimension operators can be written as a total factor Co{E/mt) for the 
leading order GF, Co{E /mt) = l — AE/Smt- The polarization function near the threshold (^) contains 
expansions in small parameters and/or v, cf. eq. (^). The leading order approximation of the low 
energy part is the exact Coulomb solution for the Green function. 

The nonrelativistic Green's function G{E) = (H — E)~^ is determined by the nonrelativistic 
Hamiltonian 

2 

H = — + V(r) (5) 
mt 

describing dynamics of the tt-pair near the threshold. The most complicated part of Hamiltonian (^) 
to find is the heavy quark static potential Vpot{f) entering into the potential V{r). The static potential 
Vpot{r) is computed in perturbation theory and can be written in the form 

Vpot{q) = ^ (6) 

that gives a definition of the effective charge related to the MS-scheme coupling constant 

av(Ai) = OLs{n){l + aiOsifi) + a2as{fif) . (7) 

Coefficients ai^2 are known [14-15]. The effective coupling is nothing but a coupling constant in 
some special subtraction scheme. The coefficient 02 allows one to find the effective /3-function (3^ for 
the evolution of the coupling at NNLO. 

High energy coefficient Ch{as) is given by the expression 

Svr \ TT J \ 27 mt J 

The NLO result for C/i(as) has been known since long ago [17-18]. At NNLO there appears a term 
proportional to the logarithm of the factorization parameter fif that separates long and short distances 
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(or large and small momenta) within the effective theory approach. The finite {^f independent) 
coefficient C2 is known ||2^, An explicit dependence of high and low energy quantities on the 
factorization scale fif is a general feature of effective theories which are valid only for a given region 
of energy. Physical quantities are factorization scale independent. In NRQCD the fif dependence 
cancels between Green's function and the high energy coefficient Ch. 

The main dynamical quantity in description of the tt system near the threshold is the nonrelativistic 
Green's function G = {H — E)~^. The Hamiltonian is represented in the form [22-24] 

H = Hc + AH, Hc = ^-^ (9) 

mt or 

with AH describing high order corrections. Constructing the Green's function is straightforward 
and can be done analytically within perturbation theory near Coulomb Green's function Gc{E) or 
numerically for complex values of E that can be used to describe the production of particles with 
nonzero width [25-33]. 

The analytical solution for Green's function is perturbative in AH 

G = Gc- GcAHGc + GcAHGcAHGc - ■■■ (10) 

Results are presented basically as an expansion in consecutive orders 

G = Go + AGi + AG2 (11) 

to check the convergence of the approximations. The leading order is given by Coulomb approximation, 
Co = Gc- At NLO the quantity AGi takes into account the corrections from the static potential 
Vpot{r) related to ai coefficient in eq. (0). At the NNLO the quantity AG2 ~ 0{al) in a sense of 
eq. (|2|) accounts for terms in the static potential Vpot{r) (02 coefficient in eq. (|7|)) and relativistic 
v'^ corrections. It also contains a second iteration of the 0{as) term from Vpot{r). The MS-scheme for 
the static potential Vpotir) is mainly used in the solution. The numerical results for GF obtained by 
different authors agree with each other [^4| . 



3 Physical results 

For the bb system an accurate description of the spectrum in terms of moments 

M^=r 

b 

is obtained in the near-threshold Coulomb PT calculations. This analysis gave the best determination 
of the numerical value for the b quark mass 

rub = 4.80 ± 0.06 GeV . 

For the tt system the top quark width Fj plays a crucial role in the calculation of the production cross 
section near the threshold. At the formal level the width is taken into account by a shift of the energy 
variable E. The mass operator of the top quark is approximated by the expression M = mt — iTt/2. 



4 



Then the kinematical variable s — 4m^ relevant to the near-threshold dynamics is substituted with 
s — 4M^ (y/s = E + 2mt) and one finds 



s - 4M^ = Amt{E + iVt) + + Ti. 

Neglecting higher orders in E and Tt one obtains a recipe for taking into account the width Ff by the 
shift E ^ E + iVt- The dispersion relation for the vacuum polarization function !!(£') has the form 



Il{E) 

With the shift recipe one finds 



p{E')dE' 
E' -E 



Because the point E + iVt lies sufficiently far from the positive semiaxis (and the origin) in the 



complex energy plane the cross section eq. (12) is calculable point-wise in energy. The hadronic cross 
section ct{E) was obtained by many authors p^]. The normalized cross section for typical values 
rrit = 175 GeV, Tt = 1.43 GeV, as{Mz) = 0.118 has the characteristic points which are usually 
considered as basic observables. They are: Ep - the position of the peak in the cross section and Hp - 
its height. In the limit of the small Tt (at least for Tt that is smaller than the spacing between the first 
two Coulomb poles) one would have Ep ~ Eq and Hp ~ 1-00(0)1^. The actual value of Tt = 1.43 GeV 
is larger than the spacing \Eq — Ei\ ~ mta'^/3 ~ 0.6 GeV therefore the peak position and height 
are not determined by the first resonance only. The convergence for Ep and Hp in the consecutive 
orders of perturbation theory near the Coulomb solution is not fast in the MS-scheme. For the typical 



numerical values of the theoretical parameters mt, Tt and as{Mz) one finds |31| 



Ep = £;o(l-h 0.58 + 0.38 -F...) 

Hp = i?o(l -0.15 + 0.12 + ...) (13) 

(see also |3^, |3^). Important contributions that affect the quality of convergence are the local term 
(~ asVo(5(r) which is related to 1/r^ non-Abelian term [^]) and higher order FT corrections to Vpot{r). 

4 Discussion and conclusion 

Slow convergence for the peak characteristics of the cross section given in eq. ( [l3[) has been actively 



discussed. The suggestions of the redefinition of the top quark mass have been made (e.g. Q^, |39|| ) 
as the use of the pole mass for a description of the cross section near the threshold is criticized on the 
ground of its infrared instability [^]. This approach is based on introduction of an effective mass that 
partly accounts for interaction [41-43]. In this talk I only discuss some possible ways of optimizing 
the convergence for the Green's function with the pole mass as a theoretical parameter Actual 
calculations near the threshold have been performed within the pole mass scheme. For optimizing 
the convergence one can use methods of exact summation of some contributions in all orders and 



renormalization scheme invariance of FT series e.g. |25]. The Hamiltonian can be written in the form 

7071^ 
4 



H = Hlo + AVpT + asVo5{f), Vq = -—^ (14) 
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where corrections are given by the perturbation theory corrections to Vpot{r) (AVpy-part) and by 
the local term (asVo(5(r)-part). The (5(r)-part is a separable potential and can be taken into account 
exactly |45]. The solution reads 

with 

Gir{E) = {Hlo + AVpT - E)-^ (16) 

being the irreducible Green's function. The PT expansion of the static potential in NRQCD is im- 
portant for getting stable results for the cross section near the threshold because the static potential 
is the genuine quantity which is computed in high order of PT in the strong coupling constant p^]. 
The convergence in the MS scheme is not fast which reflects the physical situation that the observ- 
ables represented by the cross section curve (for instance, Ep and Hp) are sensitive to different scales. 
The finite-order perturbation theory expansion of the static potential cannot handle several distinct 
scales with the same accuracy. Indeed, the PT expansion of the static potential is done near some 
(arbitrary) scale (or distance) which can be considered simply as a normalization point. The farther 
a given point lies from this normalization point the worse the precision of the PT expansion for the 
static potential is. The PT expressions in the MS scheme are not directly sensitive to physical scales 
because subtractions are made in a mass independent way (for instance, massive particles do not 
decouple automatically in the MS scheme e.g. |47, 48|). One can rewrite the static potential through 
some physical parameters which is similar to the use of the momentum subtraction scheme instead of 
the MS-scheme 

4qo 
3r 

Here tq and c are the parameters of the renormalization scheme freedom in NNLO and qq is the 
corresponding coupling in the {ro, c}-scheme [^]. They parameterize the center of the expansion (a 
normalization point) and the derivative (respective /3-function) of the static potential. The parameters 
(ro, c) can be chosen such in order to minimize the higher order corrections to a particular observable 
(e.g. ||25| where NLO analysis has been done). In such a case rg can be understood as a typical distance 
to which a chosen observable is sensitive. Note that the best approximation of the static potential 
Vpoti^) for different scales would be provided by the use of the running coupling constant as{r). The 
analytical calculation of GF becomes technically impossible in this case. The numerical calculation of 
GF requires some regularization at large distances where the IR singularity (Landau pole) can occur in 
as{r). The singularity can be dealt with if an IR fixed point appears in the evolution for the effective 
coupling constant (e.g. ||5^) or with some nonPT regularization for the potential (e.g. |^l|). For the 
top quark production the contribution of the large r region into the cross section is small because of 
the large decay width of the top quark. In the finite-order PT analysis the parameters ro and c can 
be chosen to minimize higher order corrections either to Ep or to Hp but not to both simultaneously 
because Ep and Hp are sensitive to different distances. One finds the difference of scales minimizing 
corrections to the first Goulomb resonance in NLO to be 

ln(rs/r^) = ^ + y • (18) 

Because of the large top quark width many states contribute into the position and height of the 
peak in the cross section. Therefore the characteristic distance estimates are not so transparent (the 



Vpotir) = (l + aoh Inr/ro + al{bl In^ r/rg + b'^ Inr/ro + c)) . (17) 
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NNLO peak position, for instance, is not exactly the ground state energy in the zero width limit). 
The relation ( p!8| ) can serve just as a basic guide. In practical analysis one can choose the particular 
numerical values for the parameters (ro,c) which stabilize either Ep or Hp. 
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